SHOCK WAVES FOR RADIATIVE 
HYPERBOLIC ELLIPTIC SYSTEMS 



CORRADO LATTANZIO, CORRADO MASCIA, AND DENIS SERRE 

Abstract. The present paper deals with the following hyperbolic-elliptic coupled sys- 
tem, modelling dynamics of a gas in presence of radiation, 

[ -Qxx + Rq + G ■ Ux ^ 0, 

where u G R", g G R and 7? > 0, G, L G R". The flux function / : R" ^ R" is smooth 
and such that V/ has n distinct real eigenvalues for any u. 

The problem of existence of admissible radiative shock wave is considered, i.e. exis- 
tence of a solution of the form {u,q){x,t) :— {U,Q){x — st), such that ([/, (5)(±oo) — 
{u±,0), and u± G R", s G R define a shock wave for the reduced hyperbolic system, 
obtained by formally putting L = 0. 

It is proved that, if U- is such that VAfc(u_) ■rfc(u_) 7^ 0, (where Xk denotes the k-ih 
eigenvalue of V/ and Vk a corresponding right eigenvector) and 

{£k{u-)-L) (G-rfc(ii_)) >0, 

then there exists a neighborhood U of u_ such that for any ii+ G W, s G R such that the 
triple {u- ,u+; s) defines a shock wave for the reduced hyperbolic system, there exists a 
(unique up to shift) admissible radiative shock wave for the complete hyperbolic-elliptic 
system. The proof is based on reducing the system case to the scalar case, hence the 
problem of existence for the scalar case with general strictly convex fiuxes is considered, 
generalizing existing results for the Burgers' flux f{u) = /2. Additionally, we are able 
to prove that the profile {U,Q) gains smoothness when the size of the shock ju+ — u_| 
is small enough, as previously proved for the Burgers' flux case. 

Finally, the general case of nonconvex fluxes is also treated, showing similar results 
of existence and regularity for the profiles. 



1. Introduction 

The dynamics of a gas in presence of radiation, due to high-temperature effects, can be 
modeled by compressible Euler equations with an additional term in the flux of energy. 
Dealing with small perturbation of a fixed equilibrium state in one space dimension, this 
leads to consider an hyperbolic-elliptic coupled system of the form 

(ut + f{u)^ + Lq:^ = ^^^^ 
\-qxx + Rq + G ■ Ur, = 0, 

where x G M, t > 0, -u G M", g G R and R > 0, G, L e are constant vectors. The flux 
function / : R" M" is assumed to be smooth and such that the reduced system 

ut + fiu), = (1.2) 

is strictly hyperbolic, i.e. V/(n) has n distinct real eigenvalues for any state u under con- 
sideration. For later use, we denote such eigenvalues with Ai('u) < • • • < A„(n) and with 
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ii{u), . . . , in{u), ri{u), . . . , rn{u) the corresponding left and right eigenvectors normalized 
so that li{u) ■ rj{u) = 5ij for any 

The system can be obtained from the complete gas dynamics equation with heat- 
flux radiative term by using a differential approximation of the integral equation for the 
radiative term. This approach has been proposed in the pioneering paper 0|. After that, 
equations are also called Hamer system for radiating gas (see also m for details 
on the derivation). Since the quantity q represents the radiative heat-flux term, considering 
it as a scalar quantity is physically meaningful (see for a complete physical description 
of the phenomenon). 

Under the hyperbolic rescaling {dt,dx) ^ {edt,edx), the system becomes 



Hence the rescaled hyperbolic-elliptic system can be rewritten as a singular perturbation 
of a system of conservation laws. In particular, for e sufficiently small, the system can be 
seen as a higher order correction of a viscous system of conservation laws with (degener- 
ate) diffusion term given by the rank-one matrix eR~'^L(^ G. In particular, this suggests 
that qualitative properties of solutions of (jLlj) should resemble analogous properties of 
viscous system of conservation law as soon as: u varies along a direction not orthogonal 
to both vectors L and G (non degeneracy); variations of u are mainly in the small fre- 
quencies regime (small perturbations). Indeed, failing of the first condition would imply 
degeneration of the diffusion term e R~^L G Uxx-, and failing of the second would give to 
the higher order term R~^{ut + f{u)x)^^ a dominating role. 




The above effects are present also in the scalar case, i.e. u G M (for a complete and 
introductory presentation see [23 )• As noted in ^Hj, in this situation, the 2x2 system 
enjoys many of the properties of scalar viscous conservation laws: L^-contraction, 
comparison principle, conservation of mass, constant solutions. Thanks to these properties, 
global existence and uniqueness of solutions can be proved (see for data in BV, ^H] 
for data in n ^ for the multidimensional case for data in n L°°). Nevertheless, 
regularization property does not hold: in [S] and in U^j (with a more detailed description) 
it is shown that there are initial data such that the corresponding solution to the Cauchy 
problems develops discontinuity in finite time. 

The loss of regularity appears also when dealing with the problem^ 

given u±, asymptotic states of an admissible shock wave solution to U.^) . 
does there exist a traveling wave solution to with same speed of the 

shock and asymptotic states {u±,0)? 
From now on, we refer to such a solution as a RADIATIVE SHOCK WAVE. In the scalar 
case and for f{s) = ^ s^, in ^Oj it has proved that the answer is affirmative, but that the 

-'^The analogous problem for the viscous regularization of a system of hyperbolic conservation laws is 
sometimes referred to as Gel' f and problem, and the first rigorous mathematical result has been proved in 
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profile of the radiative shock wave is discontinuous whenever the hyperbohc shock is large 
enough, i.e. — u+l is large. The precise statement will be recalled later on. In the same 
article, stability and decay rate of perturbations are determined. The absence/presence of 
jumps in small/large radiative shock wave is again a manifestation of regularity properties 
for solution of small transitions can be obtained through smooth solution, large 

transitions cannot. 

For the scalar case, many other results are available in literature. For the sake of 
completeness, let us mention them, collected in two different groups: 

- Large-time behavior, to prove stability, possibly with decay rate of the perturbations, of 
constant states jHll^, shock profiles rarefaction waves ^2]) and to find asymptotic 
profiles for such perturbations [2] ; 

- Weak solutions: to find evolution/regularity of discontinuity curves [IS], to determine 
relaxation limit under hyperbolic/parabolic rescaling \R I13j. 

The theory for system is still at the very beginning and just few results are available. 
The first paper in this direction is |7j, where global existence, asymptotic behavior and 
decay rate are proved for the solution to the Cauchy problem with initial data that are 
small perturbations of constant states. Generalizations have been given in [3] and [S], 
especially in the precise description of asymptotic profiles (diffusion waves). In particular, 
it has been proved that, for large time, the solutions to are well-approximated by the 
solutions to the viscous system of conservation laws obtained from ()1.3j) when disregarding 
the O(e^) term. The singular limit counter part, i.e. e ^ 0"*", has been dealt with in 
A different approach for analyzing the singular limit, based on the notion of positively 
invariant domain, has been considered in [22] . 

The present paper deals with the problem of proving existence of radiative shock waves 
in the case of general systems of the form p.lj) . First of all let us recall the definition of 
shock wave and radiative shock wave, where, for the sake of clarity, we refer here to the 
genuinely nonlinear case. 

Definition 1.1. A shock wave of the hyperbolic system 1)1. 2[) is a weak solution of the 
form 

u{x, t) := n_X(_^,,g) {x - st) + n+x,,^,^^) {x - st), xq^R 

where u± G M", s G M satisfy, for some k G {l,...,n}, the conditions Afc_i(n_) < 
s < Afc(n_) and Afc(n+) < s < Afc+i(n+), where Ai(n) < ••• < A„(ti) denote the (real) 
eigenvalues of V/(u) and Xiix) is the characteristic function of the set /. 

In the general case, the entropy condition for Lax shocks, Afc_i(ti_) < s < Afc(n_) and 
Afc(n-(_) < s < Afc+i(n+), should be replaced by the Liu-E condition (see the end of Section 
13 for details). 

Definition 1.2. A radiative shock wave of the hyperboUc-elliptic system (|1.1|) is a 
weak solution {u,q){x,t) := {U,Q){x — st) such that 

hrn {U,Q){O = iu±,0), 
where u± G M", s G M defines a shock wave for the reduced hyperbolic system (|1.2|) . 
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The usual Rankine-Hugoniot condition 



(1.4) 



relating the states u± and the speed of propagation s, holds also for radiative shock waves. 
This is readily seen by integrating over ah M the first equation (|1.1|) and taking in account 
the asymptotic limits of the wave. 

Being a weak solution, a radiative shock wave may be discontinuous at some .^o- Let 
{U,Q) be piecewise in a neighborhood (,^_,^+) of that is 



If the profile ?7(Co— ) 7^ f^(Co+)) then U' has a delta term concentrated at and, as a 
consequence of the second equation in the same holds for Q" and, therefore, Q is 

continuous at ^o- Hence, the first equation of suggests to call the discontinuity at 
^0 admissible if and only if the triple (U{(,oziz),s) is a shock wave for the reduced system 



Definition 1.3. An admissible radiative shock wave of the hyperbohc-elhptic system 
is a radiative shock wave (U, Q) if there exists a discrete set J C M for which U is 
off^ J', Q is continuous on M and off and the triple {U{(,o±), s) is a shock wave 
for dOJ) for any G J- 

Later on, it will be shown that the admissible radiative shock wave we deal with have 
at most one discontinuity point. 

Now we are in position to state the main problem we deal with: 

Problem. Given a triple u± € M", s G M, defining a shock wave for the 
reduced system i)1.2|) . does there exist a corresponding admissible radiative 
shock wave for the hyperbolic-elliptic system i)J. J|) ? 

The analogous problem has been affirmatively solved for other regularizations of (|1.2|) 
(see |17| and |^ for the viscous and the relaxation approximation, respectively). 

In the context of radiative gas model, this problem has been addressed for the scalar 
case ^lE] and for specific systems [T^ . 

The result in jl9j concerns the existence of radiative shocks for convex fluxes /. Ele- 
mentary computations shows that the parameters e and m in the cited paper are related 
with L, G, R by the relations e = LG/R and m = ^/R/LG. Hence the existence result in 



QeC\{C^,Co])nC\[Co,^+)). 



can be stated as follows. 
Theorem 1.4. [1^] Assume /" > 0. // 




(1.5) 



then there exist a continuous radiative shock. 



^Given J" C R, we say that a function F : R ~+ R is C''' OFF J if for F G C'^{[a,b]) for any interval 
[a, b] such that J n (a, b) is empty. 
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Condition (|1.5j) is satisfied in the case of small shocks, i.e. if |m_ — "U+I is small. Also, 
by a hyperbolic rescaling ^ L~^dx, dt i— > L~^dt, the system becomes 

ut + f{u)x + qx = 

where m} = R/Lp'G^ and v = GL/R. As m ^ 0, formally, we get a scalar viscous 
conservation law. Hence condition (|1.5j) can be read as a measure of the strength of the 
viscosity term, needed for smoothness. 

In [TUl, it has been considered the scalar case with the choice of a Burgers' like flux 
function f{u) = ^ and, without loss of generality, L = R = G = 1. 

Theorem 1.5. ^Hl Let u± be such that n+ < u_ and set s := {u-^ + n_)/2. Then there 
exists a (unique up to shift) admissible radiative shock wave {U,Q). 

The authors can also determine the smoothness of the profiles {U, Q) as the size |u+ —u^ \ 
varies, showing that regularity improves as the amplitude of the shock decreases: 
-if |u+-u_| > ^^2, then U E C°(M\{Co}) for some ^0^1^ and Q is Lipschitz continuous; 
- if \u+ - M_| < 2V2^/(n + 1) for some n G N, then U G C"(M) and Q G C"+i(R). 

In the case of systems (i.e. u G M"), the literature on the subject restricts to the very 
recent article ^SJ- There, the authors consider the problem of existence of smooth radiative 
shocks for a specific model describing the gases not in thermodynamical equilibrium with 
radiations. In their case, the coupled elliptic equation is nonlinear with respect to the 
temperature. 

In the present article, we deal with general systems and consider admissible ra- 
diative shock waves, hence possibly discontinuous. We are able to prove that, for small 
amplitude shock waves of H1.2() . there always exists a radiative shock profiles. 

Theorem 1.6. Let U- G M" be such that the k—th characteristic field of H.<^) is genuinely 
nonlinear at n_, that is VAfc(u_) • r}^{u^) ^ 0. Assume that 

{£k{u^)-L){G-rk{u^))>Q. (1.6) 

Then there exists a sufficiently small neighborhood U of ti_ such that for any uj^ G li, 
s G M such that the triple (ii_,u+; s) defines a shock wave for there exists a (unique 

up to shift) admissible radiative shock wave for 

The proof is essentially divided in two steps. The first one is to reduce the system 
case to the scalar case. This (surprising!) possibility is essentially due to the fact that 
the diffusion matrix L ® G \s rank-one. The second step is to generalize the existing 
proof for Burgers' like flux to general strictly convex fluxes and to show the existence 
of an heteroclinic orbit connecting the asymptotic states. The analysis is complicated 
by the fact that the differential equation for the proflle is not in normal form and that 
discontinuity of U can arise. In dealing with this problem, it turns to be useful to work 
with the integrated variable z such that q = —Zx, hence to deal with the equations 

ut + f{u)x = LZxx ^^^^ 



6 



CORRADO LATTANZIO, CORRADO MASCIA, AND DENIS SERRE 



instead of the original system 

Remark 1.7. The smallness assumption on the shock u+; s) is needed only in the first 
step of the proof, that is in the reduction procedure from system to scalar case, which 
is carried out in terms of Implicit Function Theorem. Since we are able to prove the 
existence of (possibly discontinuous) radiative shocks for general strictly convex scalar 
models and general large admissible shocks, our existence result includes in principle the 
case of discontinuous radiative shocks for 

By applying a similar strategy to the one used in ^0], we can also get analogous result 
on smoothness of the profile (U, Q). 

Theorem 1.8. There exists a sequence {Sn}, ^ as n ^ oo, such that the profile u 
is C^+^ whenever — < £«. 

For completeness, once the convex case has been treated in full generality, we consider 
the case of nonconvex fluxes, again showing the existence and regularity of admissible 
radiative shock profiles. The number of possible jumps of the profiles and the construction 
of the profiles themselves are strictly related with the number of inflection points of the 
flux function (see Section 0)) . 

The paper is organized as follows. Section |21 is devoted to prove the existence result in 
the scalar case for strictly convex flux functions /. Gaining of regularity of the profile as 
the amplitude decrease is considered in Section |S1 again in the strictly convex case. The 
scalar case with a general (nonconvex) flux function / is treated in Section^ Finally, in 
Section |S1 we show how to reduce the problem from the system case to the scalar one. 



We start our analysis with the case of a scalar conservation law, with convex flux, 
coupled with the linear elliptic equation describing the radiating effects. Specifically, we 
discuss the existence of a travelling wave profile for the 2x2 hyperbolic-elliptic system 



where x E M, t > 0, n and q are scalar functions, the flux function / is strictly convex, 
R > and L, G are constants such that LG > 0. As reported in the Introduction, previous 
results on the same problem are contained in (general convex fluxes, smooth radiative 
shocks) and in jJU] (Burgers' flux, general radiative shocks). 
Applying the rescaling 



2. Existence of the profile in the convex case 





(2.1) 
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We stress that, thanks to the discussion of Sectional the results of the present section 
win give the existence of a /c-travelhng wave solutions for system with genuinely 

nonlinear k characteristic field. The study of the general scalar case is left to Section 0] 
and the results for system Hl.lj) without GNL assumptions is again a consequence of the 
reduction arguments of Section [51 

Remark 2.1. The existence of the profile does not need the strict convexity of the flux 
function /. Indeed, our proof applies also to the case of a function / such that f{u) — cu 
is monotone on two intervals and (n^,,n_). 

Let us consider a solution to H2.1|) of the form {u, q) = {u{x — st), q{x — st)) such that 

/ , N fiu+) - f{u-) 

u(±oo) = u±, s = , n_|_ < 

ii+ — u_ 

We introduce the variable z as the opposite of the antiderivative of q, that is z^ := —q, 
with z(iboo) = z± = u±. Hence, after an integration of the second equation, system ()2.1|) 
rewrites as 

\ut + f{u)x - z^x = 
\-Zxx + z-u = {). 

Thus, the equations for the profiles u{x — st) and z{x — st) take the form 

j-su' + f{u)' - z" = 
\-z" + z-u = Q, 

that is, after integration of the first equation 

f -s{u - u±) + /(u) - f{u±) = z' ^2 2^ 

]u = z — z" . 

At this point, we rewrite H2.2|) as the following second order equation for z 

z' = F{z - z";s), (2.3) 

where F{u; s) = f{u) — f{u±) — s{u — u±). We shall prove the existence of a profile for 
z solution of (|2.3j) between the states z_ = u_ and z^ = u^^ z^ > z+, which will give 
the existence of our profile for u = z — z" . To this end, let us note that, thanks to the 
strict convexity of /, the function F(-;s) is strictly decreasing in an interval [2;+,2;*], and 
strictly increasing in where F{z±;s) = and F(z*;s) = —m < 0. Hence, F{-;s) 

is invertible in the aforementioned intervals and we denote with h± the corresponding 
inverse functions. In the next proposition, we analyze the behavior of the two ordinary 
differential equations z" = z — hj-{z') in the corresponding intervals of existence. 

Proposition 2.2. Let us denote with zj^ = zj^{x) the (unique up to a shift) maximal 
solution of 

z" = z-h+iz'), (2.4) 

with z^{+oo) = z^ and z'_^{+oo) = 0. Then z^ is monotone decreasing, z'_^ is monotone 
increasing and moreover Zj^ is not globally defined, that is, there exists a point, assumed 
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to be (thanks to translation invariance) ,such that 

z+(0) - 4(0) = z*, z^(0) = -m. (2.5) 
Proof. We rewrite (|2.4|) as a system of first order equations as follows 

X' = H+[X), (2.6) 

where 



Then 



X-( J, i7+(X)- , 







and (^;+,0) is a saddle point. We are interested in the stable manifold at this point and, 
since we want z to be decreasing, let us follow the trajectory which exits from (z+,0) in 
the lower half plane of the states space. Then we claim that z is monotone decreasing and 
z' is monotone increasing. Indeed, if by contradiction z attains a local maximum, say at 
X = xq, then z'{xq) = and > = z — /i+(0) = 2 — z+ at x = xq, which is impossible. 
Thus z is monotone decreasing and z' < 0. Now, if by contradiction z' attains a local 
minimum at x = xq, then the trajectory z' = (p{z) in the {z,z') plane must attain a local 
minimum at x = xq, that is <f'{z) = and <f"{z) > 0. Thus, at x = xq, 

-h+iz') 



= ^'{z) = - 



z' 



and 



, , , d [ z-h+{z') \ {z'-h\{z')z")z'-{z-h^{z'))z" 1 
^'^ = d-z{ 1' )= = I^<°' 

which is impossible. Hence z' is monotone increasing and clearly hj^{z') G [z_|_,2,,]: There- 
fore z" = z + 0(1), which implies the solution does not blow up in finite time. 
Thus, we are left with the two following possibilities: 

(i) the solution is defined for any x G M; 

(ii) the solution reach the boundary of the domain of definition of the differential 
equation in a finite point, which can be assumed to be up to a space translation, 
namely, 1)2. 5|) holds. 

Hence, the proof is complete if we can exclude the case (i). To this end, let us assume that 
(i) holds. Thus in particular z' is bounded for any x S M, because z' G {—m, 0). Since z is 
monotone decreasing, then z ^ I £ (z+, +00] as x — > — cxd. If / < +00, since z' is increasing, 
then z' — > as X ^ —00 and therefore z" = z — h^{z') ^ / — /i+(0) = I — z^ = 0, which 
is impossible. On the other hand, let us assume / = +00. Since z' is globally bounded, z 
decreases at most linearly, which is incompatible with z" = z + 0(1), because the latter 
implies an exponential rate. Therefore, the solution cannot be defined for any x G M and 
the proof is complete. □ 
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With the same kind of arguments, it is possible to analyze the behavior of the maximal 
solution z_{x) of the equation z" = z — h_{z'). Hence, the following proposition holds. 

Proposition 2.3. Let us denote with Z-{x) the (unique, up to a space shift) maximal 
solution of 

z" = z-h^{z'), (2.7) 

with z^{—oo) = Z- and z'_{—oo) = 0. Then z^ and z'_ are monotone decreasing and 
moreover Z- is not globally defined, that is, there exists a point, which we can assume to 
be (thanks to translation invariance), such that 

z_{0) - zl{0) = z^, z'_{0) = -m. 

With the aid of Propositions 12.21 and 12.31 we shall build a trajectory joining z^ and 
z-t- by finding a point of intersection for the orbits of the maximal solutions z-{x) and 
z+{x) in the state space {z,z'). The monotonicity of z±{x) and z'j.{x) shall guarantee 
such an intersection is unique and therefore the resulting trajectory from z- to is 
unique, up to a space translation. The existence of the aforementioned intersection is a 
straightforward consequence of the following lemma. 

Lemma 2.4. Let us denote with z+{x) and z^{x) the maximal solutions of \2.4\l and \2. 1\) 
respectively. Then 

2+(0) > > z-(0). (2.8) 

Proof. We prove only the inequality on the left of ()2.8() . the right one to be proved similarly. 
Let us denote with y the heteroclinic orbit of 

{y' = F{y-s) 

\y(±oo) = z±. 

Our aim is to compare this solution with the solution Zj^{x) in the phase space. To this 
end, as before we denote with = ip{z) the function whose graph is the trajectory of 
{zj^{x) , z'_^{x)) in the {z,z') plane and with tp = ip{z) = F{z; s) the function whose graph 
is the trajectory of {y{x),y'{x)) in the {z,z') plane. Then 

iP'iz+) = F'{z+-s)<0 

and V''(^+) = '^i is the negative eigenvalue of V-ff+(z+,0), with defined in the proof 
of Proposition 12.21 namely, the negative root of 

Since P(F'(z+;s)) = (F'(2;+;s))^ > 0, then F'{z^;s) < Ai and therefore the trajectory 
iIj{z) leaves the z-axis in below the trajectory f{z). 

If ip{z) intersect ip{z) in a point of the (z, z') plane, then in that point we have (p{z) = 
^{z) and ip'{z) < Tp'{z). Hence, there exist a point x such that y{x) = -z+(x) and y'{x) = 
z',(x) < 0. Moreover, since 
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and 




in that point we also have y" < z"^ = zj^ — h^{z'_^) = y — h^{y') = 0. On the other hand, 
> y" = F'{y;s)y' = F' (y; s)F{y; s), that is F'{y;s) > because F{y;s) < 0, namely 
y = 2+ > in that point. 

Conversely, let us assume that (p{z) remains above ip{z). From Proposition 12 . 21 we know 
that z'_^{0) = — m, which is the minimum of ip{z), attained at z*. Therefore we must have 
2;+(0) = z^ and the proof is complete. □ 

With the aid of Lemma 12.41 we are able to prove the main result of this section, which 
is contained in the following theorem. 

Theorem 2.5. There exists a (unique up to space translations) profile z with z{±oo) = 
z± = u± such that the function z{x — st) is solution of i2.!A) . where s is given by the 
Rankine-Hugoniot condition. The solution z is away from a single point, where z" has 
at most a jump discontinuity. 

Moreover, there exist a (unique up to space translations) profile u with u{ztoo) = u± 
such that the function u{x — st) is solution of h2.^3\) . where s is given by the Rankine- 
Hugoniot condition. This profile is continuous away from a single point, where it has 
at most a jump discontinuity which verifies the Rankine-Hugoniot and the admissibility 
conditions of the scalar conservation law ut + f{u)x = 0. 

Proof. Let us observe that Lemma 12.41 implies that there exists a point in the {z,z') 
plane where the graphs of z+ and z- intersects, namely z_ = 2;+ and z'_ = z'_^ in that 
point. Moreover, due to the monotonicity of these graphs, which comes from the proofs 
of Proposition 12.21 and Proposition 12.31 this intersection is indeed unique. Hence with 
appropriate space translations, we can find a point x such that Z-{x) = z+{x) = z and 
z'_(x) = z'j^{x) = z. Set 



Thus, this profile is the unique (up to space translations) solution of ()2.3() and it 
verifies z(iboo) = z±. Moreover, z{x) is in the intervals (— cxd,x] and [3;,+oo) and 
finally z"{x - 0) = z'L{x) = - h^{z'_{x)) = z- h-{z) and z"{x + 0) = ^(x) = 

z+{x) — h+{z'_^_{x)) = z — h^{z). 

The regularity of u = u{x — st) is a direct consequence of the first part of the theorem 
and of the relation u = z — z". Moreover, in the case of a discontinuity in n, namely 
u{x — 0) / n(x + 0) , {u{x — 0) , u{x + 0); s) verifies the Rankine-Hugoniot condition for the 
strictly convex conservation law ut + f{u)x = 0. Indeed, u{x — 0) = h^{z), u{x+0) = h^{z) 
and a direct calculation shows 



Finally, u{x — 0) = h^{z) > h^{z) = u{x + 0), that is, this shock is admissible and the 




z^{x) X < X 
z+{x) X > X. 



f{h+{z))-f{h^{~z)) 
h+{z) — h^{z) 



= s. 



proof is complete. 



□ 
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Remark 2.6. The above theorem contains an uniqueness result, among the class of radia- 
tive shocks given in Definition 11.21 An uniqueness result in a wider class of solutions is 
contained in Theorem 14.31 We postpone this result because it is valid for general fiux 
functions, disregarding convexity properties. 

Remark 2.7. From the last lines of the above proof, it is clear that the profile u can have a 
jump discontinuity in a point, as it was already proved in ^U] for sufficiently large shocks 
in the case of the Burgers' equation. Moreover, in that paper it is proved also that, below 
an explicit threshold, the profile is continuous and it smoothes out as the strength of the 
shock decreases. In our case, the profile in u is continuous if z"{x — 0) = z"{x + 0), namely 
when h-{z) = = z*, that is —m = z = z'j^{x), z±{x) = and z'j^{x) = 0. This 

property and the further regularity of the profile for sufficiently small shocks is proved in 
Section |21 below. 



3. Regularity of the profile in the convex case 

In Section|21 we proved the existence of the travelling wave (n, q) = {u{x — st), q{x — st)), 

(, N f{u+)-f{u-) 
u[±oo) = u±, s = , 

li+ — U- 

for ()2.1() . when the fiux f{u) is smooth and strictly convex; in the present section we focus 
our attention to its regularity. 

Let us start by observing that this issue is related to the smallness of the shock we are 
dealing with jlUj . Here we shall prove that it is possible to recast such a property in terms 
of smallness of the term —qxx in (|2.1j) . namely when this system increases its diffusive 
nature, being closer to its parabolic equilibrium. This will be made in terms of a diffusive 
scaling with respect to the shock strength e = \u^—u- \ and, following the ideas of ^Hl) by 
analyzing the system for (n, v = u'). In that way, we will obtain the same kind of results of 
jlUj for the special case f{u) = ^u^. Let us finally observe that similar phenomena arise in 
the discussions of Section to prove the existence and regularity of travelling wave profiles 
of (|4.H) with fiux functions / with change of convexity. In that case, the small parameter 
e is already present in the model and it is not connected with the smallness of the shock. 
However, existence and regularity of the profile will require once again the smallness of 
that parameter, namely, as before, when the hyperbolic-elliptic model is close enough to 
its diffusive underlying dynamic. 

From 1)2. 1() and after integration with respect to ^ of the first line, the equations for the 
profiles u and q are given by 

-s(n - u±) + f{u) - f{u±) + q = 
u' = q" - q, 

because q{±oo) = 0. Thus, using q = — [f{u) — f{u±) — s{u — u±)] 
we end up with 



(3.1) 

= -F{u;s) in 



u' = -F{u; s)" + F{u- s) = -{f'{u) - s)u" - f"{u){u'f + F{u; s). 
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Hence, we obtain the following system in the state space (w, v = u') 

(3.2) 



u' = V 



[{f'{u) - s)v' = -f"{uy -v + F{u- s). 

It is worth to observe that system (|3.'2() is singular where {f'{u) — s) = and, thanks to 
the strict convexity of /, the latter occurs at a unique u G At this stage, driven 

by ^H], we scale both the independent and the dependent variable as follows: 

Then, dropping the tildas, system (|3.2|1 becomes 

{u' = V 
f'e{n)v' = ^ i-e'f^'iuy - V + Mu)) , ^^'^^ 

where fe{u) := -^F{uj^ +£U]s) and u G [0,1]. Once again, system (|3.3|) is singular at 
u = vf , for a unique G (0, 1). Finally, as in JO]) we remove the singularity of (|3.3|) by 
using once again a new independent variable r] defined by 



oo 



V 



The resulting system reads as follows 
in' = n{u)v 

\v' = ^{-e^n'{uy-v + Mu)). ^ ■ ' 

We notice that (|3.4|) admits (0,0) and (1,0) as equilibrium points, corresponding to the 
two equilibrium points (u±,0) of the original system. 

Remark 3.1. (i) For any n > 2, if the original flux function / is C", then — > i/"(u+)u(u— 
1) in C"([0,1]), asejO. 

(ii) The last change of independent variable gives a reparametrization of the orbit of ()3.3|) 
in the two regions [0, n^) and {u^, 1]. Let us first observe that f^{u) < for < n < 
and f^{u) > for < u < 1. Hence, we obtain continuity of the orbit u of ()3.3() . provided 
we prove the two previous reparametrized orbits verify u{—oo) = 0, ti(+oo) = u^, and 
u{+oo) = u^, u{—oo) = 1 respectively. We shall prove the former property, the latter 
being similar. 

(iii) As for 10 , further regularity for the profiles solutions of l|3.4j) implies further regular- 
ity of the original profile, thanks to their exponential decay toward the asymptotic states 
at iboo. 

For the sake of clarity, we state first the result concerning the continuity and the 
regularity for the profile u solution of 1)3. 4() and then the one for the C"^ regularity. The 
general case, including the former ones, is then stated in the last Proposition. Here below, 
we shall assume / to be smooth. 
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Proposition 3.2. There exist two values e, > such that, fore < min{e, eo}; the orbits 
of which pass through the equilibrium points (0,0) and (1,0) meet at the equilibrium 
point {u^,V2), where 



_ -l + v/l + 4e^/i-(^^)/,(n^) 
In particular, the orbit u is . 

Proof. We start by studying the "new" equilibrium points of ()3.4|) introduced by the last 
change of variable, besides the aforementioned points (0,0) and (1,0). Hence we must 
satisfy the relations 

f/i(«) = o 

\e''f':{u)v^ + v-fe{u) = (), 

that \s u = and 

_ -l^^Jl-Ae^f'^{u-)\f,{u-)\ 

provided £ < £o, thanks to Remark l3.1H (il. Moreover, using again Remark l3.H -fi). there 
exists an e > such that for any e < e the whole curves 



_ -lTVl-4e^f^'{umn)\ 
^i,2(^)- 2e2/^'(n) ^"^-^^ 

with u G [0, 1] are real. We shall now study the nature of these points to prove the first 
assertion of the theorem. The Jacobian associated to system (|H.4j) is given by 

/ f^\u)v f'^u) 

J{u,v) = 

\-n"iu) + ^/i(n) -2f^'{u)v - - 

Evaluating J(0,0), J(1,0) and J{u^,V2) we conclude (0,0) and (1,0) are saddle points, 
while (ti^,^!) is a sink. Moreover, the curve v = V2{u) given in (|ci.6|) passes through the 
point (0, 0) and 



= m- 

u=0 

Let us consider the orbit exiting from (0,0) at —oo, namely the one on the unstable 
manifold of that (saddle) point, in the halfplane v < 0. Its tangent vector is then given 
by the eigenvector related to the positive eigenvalue A+ of J(0,0), that is 




\fmj 

A direct calculation shows 

namely the trajectory leaves the equilibrium point above the curve v — 
a straightforward analysis of the dynamical system (|3.4j) implies it must reach the equi- 
librium point {u^,V2) at +cxd. The same result holds for the orbit exiting from (1,0) 



14 CORRADO LATTANZIO, CORRADO MASCIA, AND DENIS SERRE 

and therefore the first claim of the theorem is proved. Finally, since v is continuous and 
bounded everywhere, u £ and the proof is complete. □ 

Proposition 3.3. There exists a value < ei < eo such that for e < mm{e,ei}, the orbit 
u of 1^ is C72. 

Proof. We evaluate the Jacobian J{u,v) at the equilibrium point (u^,?;!) to obtain 
/ f^'mvl 

J(n^^;i) = 

\-f^"{u^)mf -2f^'(unq-- 

which admits the following (negative) eigenvalues 



■^1 — Je 1^2 — 



For e = eq, A| = and Af = < 0, while A| -oo and Af -f"{u'^)\fe{u'^)\ < as 
e I 0. Thus, there exists an ei < Eq such that for e < ei, Af < Af < 0. Therefore, for any 
e < ei, the orbit in (u^,?;!) is tangent to the eigenvector related to Af, that is 

1 




instead of the eigenvector related to A|, namely (0, 1)* 
Finally, 

dv 



du 



^ ,3.7) 

which implies the regularity of the u-component of the orbit of 1)3. 4|) . that is, the 
regularity of the n-component of that orbit. □ 

To prove further regularity of the profile {u,v) solution of (|3.4|) . we follow again the 
ideas of ^U]. For n > 0, let us consider by induction the new variable 



Wn = 7 I V 



n 



Y,y^'^{n-uy I , (3.8) 



that is 



u = ^ w^j{u - uy + Wn{u - n^)" := /i„(n, Wn). (3.9) 

j=0 

Let us note that, with the notation introduced above, we choose Wq = wf given by 

p.7|) . Hence, for n = 1, 1)3. 8|) reduces to 

V -V2 _e 
Wl = — - W^, 

u — u'^ 

that is, u)f is the value of the first derivative of the curve v = v{u) for u = and therefore 
the rest point {vF ,V2) becomes (n'',0) in the new pair of variables {u,wi). However, for 
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n > 1, we assume by induction the constants Wj, j = 1, . . . ,n — 1 in (|H.8|) and (I.S.tjj) are 
given, but we determine by imposing a different condition (see (|3.14|1 below) and we 
obtain the relation 
1 



for any j >2 



as a consequence. The system for (u, Wn) is given by 

^u' = f'^{u)hn{u,Wn) 



(3.10) 



w. 



e2(ii - u^y^ 



(3.11) 



where 



{u, Wn) = -e'^f'J{u)hn{u, Wnf - hn{u, W^) + fe{u) 



We are ready now to prove the main result of this section, namely, Proposition l3.4[ (already 
proved for n = 0, 1 in Propositions 13.21 and 13. 3|) . 

Proposition 3.4. For any n > 0, There exists a decreasing sequence of positive values 
{e^n}n>0; such that, for any e < min{e, e„}, v is a C" function of u and admits the 
expansion 

n 

v = Y^ w'j{u - u'y + o{{u - -u^)"), (3.12) 

j=0 

for u ^ . 

Proof. A direct calculation shows 

9n{u,Wn) = Fn{u) + G„{u)Wn + Hniu)wl, 

where 



Y,^',iu-uy+feiu) 

j=0 



.j=l I \/=0 



Gn{u) = -2e^f'jiu) E^K^ - ^y^^ I - - ^"y 

\3=0 

-e^feiu) \^{3+n)w]{u-u'y+^-^ | - e^nf',{u)wl{u - u'^^-^- 
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Hn{u) = -e^ {fj{u){u - u'f^ + nf,{u){u - n^)^-!) . 
In addition, for sufficiently smooth fluxes /, since f'^{u^) = 0, we obtain, for u 
Gn{u) = 0{{u - n^D, Hn{u) = 0{{u - u'f^). 



Moreover, the expression of Fn{u) implies Fn{u) = F„_i(n) + 0{{u 
the coefficients of the Taylor approximation of Fn{u) about if do not depend on n, that 
is 

n 

F^{u) = cjiu - uy + 0{{u - uT^^), 

j=0 

with Cj independent from n, for any j. We write down these coefficients more explicitly 
as follows 

Cj = aj{wQ, . . . ,iDj_^) + /3j{tDo)tDj, 

where 

Pj{wl) = -{e\2 + j)f'Jiu^)wf, + l). 
Now, using again Remark l3.H -(i). there exists a decreasing sequence {ej}j>o such that 

/3,K)<0 (3.13) 
for any e < ej, and therefore we choose iUj, for e < ej and for any j, such that 

cj = ^ F^i\u^) = 0. (3.14) 

It is worth to observe that, if we characterize the above argument, and in particular the 
values Sj and 1)3. 14() . to the cases n = and n = 1, we recover the previous choices Wq = 
and w)f given by (|3.7() . 

At this point, we have determined all coefficients tD^, j > in (|3.12|) and we are left to 
the proof of the regularity of u as a function of u close to n = -u^. Assume by induction 
the result holds for n and introduce the new variable Wn defined in 1)3. 8() . with Wj, j < n 
verifying (|3.14|) . Therefore, Wn as u ^ hy the induction hypothesis and, thanks to 
(|3.14() . (^^,0) is a rest point for ()3.11() . which correspond to the original rest point (■u'^,^!) 
for H3.4|) . Since hn{u^ , Wn) = for any Wn and in view of the structure of 9niu, Wn) showed 
before, the Jacobian of 1)3. 11() evaluated at this rest point is given by 

/ r:{n')wf, \ 



^. Fl^{u){u - u') - nFn{u) Gn{u) 
lim — ^ — lim 



(3.15) 



Since we know (|3.14|) for j < n, then 

u^u^ e^(ti — ti^)"+^ 
Moreover 
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Therefore, the eigenvalues of (|l'{.15j) are given by f"{u^)wQ = Af and Al — nAf with eigen- 
vectors 



and 



V(n + l)Af-Aiy 

Thus, the result for n + 1 is obtained proceeding as in Proposition 13.31 by imposing 
Af > Al — nXi, that is, condition (|3.13|) for j = n + 1. Finally, from the uniqueness of the 
Taylor expansion of v about u^, relations (|3.1fl|) is also verified. □ 

Remark 3.5. In the case of a Burgers' flow, that is, the case studied in jlOj . our scaling 
procedure will lead to a modified quadratic flux independent from e, namely ^u{u— 1) (see 
also for further comments on the behavior of the flux with respect to this scaling). 
Therefore, for that particular case, the thresholds £j can be explicitly calculated from 
1)3. 13p , e = eo and our results coincide with the ones established in [TUl . 

4. The non convex case 

Now, let us consider existence and regularity of travelling wave profiles when the flux 
function may change its convexity. As already pointed out in Section |31 in that case, we 
shall prove these results if the behavior of ()2.1|) is close enough to its underlying diffusive 
equilibrium, that is, when —Qxx is sufficiently small. Hence, let us consider the 2x2 system 

iut + f{u)x + qx = ^^^^ 

with < e ^ 1. Let (n_,u+;s) be an admissible shock for the inviscid conservation law 
ut + f{u)x = 0, namely, it satisfies the Rankine-Hugoniot condition 

f{u+)-f{u.) 
s = 

U+ — U- 

and the the strict Oleinik condition 

fju) - f{u^) ^ f{u+) - f{u^) ^ f{u+)-f{u) ^ 

u — U- M+ — n_ Uj^ — u ' 

for any u between n_ and u+. As before, let us introduce a new variable z as follows 
—Zx'-=qso that the equations for u{x — st) and z{x — st) are given by 

|^(";;1=^' (4.3) 

I —£z + z = u, 

where, as usual, F{u\ s) = f{u) — f{u±) — s{u — u±) and z{±oo) = z± = u±. Without loss 
of generality, assume n_ > n+. Then (|4.2[) becomes 

F{u;s) < 0, for any u £ {u+,U-) (4.4) 

and in addition it implies the Lax condition 

/'(n+) < s < f{u^). 



18 



CORRADO LATTANZIO, CORRADO MASCIA, AND DENIS SERRE 



We start by considering the non degenerate situation, that is 
/'(n+) < s < f'{u^), 

, (4-5) 
/ {u) 7^ for any u G (n_|-, n_) with / {u) = 0. 

The result in the general case will be proved by an approximation procedure at the end of 
the section [2111 • shall construct a profile for the z-component by solving the equation 



z' = F{z - ez"- s) (4.6) 

and then u will be given hy u = z — ez" . 

Since / is not strictly convex, in order to invert the function F(-;s), we decompose 
in the disjoint (up to vertices) union of subintervals /i,...,l2n, n > 1, where 
F{-;s) is monotone. More precisely, thanks to 1)4. 4|) and ()4.5() . F{-;s) is decreasing in 
l2k-i and increasing in l2k, k = 1, . . . ,n. In the spirit of Propositions 12.21 and [731 we 
shall construct maximal solutions zq, z\, zJ,, k = l,...,n — 1 and z„ that correspond 
respectively to /i, l2k, hk+i^ k = 1, . . . ,n — 1 and /2n and then match the corresponding 
graphs in the phase plane to obtain a global profile as in Theorem 12.51 In addition, 
this solution will be in the points {z2fc} = ^2k H /2fc+i7 k = 1, . . . ,n — 1, where F{-; s) 
attains a local maximum, while the second derivative will have a jump discontinuity in 
the points {^^2^-1} = hk-i H hk, k = 1, . . . ,n, where F{-; s) attains a local minimum (z|, 
k = 1, . . . , 2n are the zeros of F'{-; s) in z-], ordered from the left to the right). In 
this way, the resulting profile for u will be regular, except for n points, where it has at 
most a jump discontinuity, which is an inviscid shock satisfying 1)4. 4() . Clearly, the the case 
n = 1 correspond to the strict convex case treated before. Moreover, the construction of 
the (unique, up to a space translation) maximal solutions z^ and z_ of Propositions 12.21 
and 12.31 can be repeated under assumptions 1)4.4(1 and (|4.5() to obtain the (unique, up to a 
space translation) profiles zq and z^ corresponding to the intervals z^] and [z2n-i' ^-]- 
Taking into account also the results of Lemma 12.41 we know that such maximal solutions 
verify, up to a space translation, the following properties: 

(i) Zq : [0, +00) (z+, 2:0(0)], Zq monotone decreasing, z'q monotone increasing and 

2:o(+cx)) = z+; 

(ii) Zn : (—00,0] [z„(0),z_), Zn and 2;2n monotone decreasing and Zn{—oo) = z^; 

(iii) 2:0(0) > zl and 2:„(0) < 22n-i- 

It is worth to observe that the existence of the above profiles with the aforementioned 
properties does not depend on the value of e > 0. Hence, we are left with the construction 
of the intermediate maximal solutions z^ i and Zk^T^ /c = l,...,n — 1, solutions respectively 
of 

ez" = z - h2k{z') (4.7) 

and 



ez" = z - h2k+i{z'), 
where hi denotes the inverse of F{-] s) on Ii, i = 1, ... , 2n. 



(4.8) 
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Proposition 4.1. Let us assume conditions i4.4\ ) <ind i4.5\ ) hold. Then, for any k = 
1, . . . ,n — 1, there exists a (unique up to space translations) maximal solution z^^i of \4-'l\j 
and Zk^r of i4-^ > with initial data 

( zk,i{o) = zkAo) = 4k (.c,) 

\4,iiO)=z',,M=Fiz;,;s). ^ ■ > 

Moreover 

4/(0) = 4,(0) = (4.10) 

and 

(i) Zk^i, z'j, I are monotone decreasing and z^^i is not globally defined, that is there exists 
a point S^k^i > such that 

Zk,i{lk,i) - ez'l,i{ik,i) = 4fc-i> 4,/(6,/) = F{z*2k-i]s)] 

(ii) Zk^r is monotone decreasing, , monotone increasing and Zk,r is not globally de- 
fined, that is there exists a point ^k,r < such that 

Zk,riCk,T) — ^Zk,r(.^k,r) = ^2fc+l' ^•,r(?fc,r) = -^(^2/c+l ! '5) ! 

(iii) there exists an such that, for any e < e^, 

ZkMk,l) < Z*2k~l, Zk,r{^k,r) > Z2k+1- (4.11) 

Proof. Let us start by justifying tiie clioice of initial data in ()4.9|) . Since we want to joint 
Zf^^i and Zk^r in to obtain a smooth (say, C^) profile, we are forced to choose H4.9|) 9. 
Moreover, differentiating (|4.6() we get 

z" = {z- ez!')'F\z - ez"- s). (4.12) 

In the sequel, we shall also need 4/(0) 4'r(0) finite and therefore, assuming that 
requirement, form 1)4. 12() we obtain (|4.10|) . which implies H4.9() i . in view of l|4.7() and (|4.8() . 

Since F{zl^;s) < 0, we can extend continuously the function h2k at the right of 
to obtain the existence of a (not necessary unique) maximal solution of ()4.7|) - H4.9() . 
which coincides with ^a:,/(?)> fo^' ^^y ^ > 0, as long as Zk^i{S,) < z^j^. Hence, property (i) is 
proved as before: /i2fc is the inverse of F{-; s) on an interval where F{-; s) is increasing and 
therefore it corresponds to the case studied in Proposition 12.31 A symmetric argument 
will lead to the situation of Proposition 12.21 which gives (ii) . The remaining part of the 
results will be proved only for Zk^i, the ones for Zk^r being similar and thus we drop the 
subscript k,l. 

Let us now prove 

\z"'{0)\ < +00. (4.13) 

To this end, let us come back to the original variables u and q. From ()4.1() . and after 
integration with respect to ^ of the first line, the equations for the profiles u and q are 
given by 

] F{u;s) = -q 



u' = eq" 
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because q{±oo) = 0. Thus, the profile u, regular away from and continuous there, verifies 



u 



' - -eF"{u; s) + F{u; s) = -eF'{u; s)u" - eF"{u; s){u'f + F{u; s) 



Hence, we obtain the following system in the state space (n, v = u') 



u' = V 



[eF'iu; s)v' = -eF"{u] s)v'^ - v + F(n; s). 
At this point, as in JOl) we shall introduce a new independent variable r] given by 

F'{u{0;s)dQ, 



(4.14) 



which will move the singularity attained for ^ = to r/ = oo, thanks to exponential decay 
of F'[u{rf)]s) toward zero, as r/ ^ cxd (the same kind of procedure used for regularity 
results in Sectional). Hence, (|4.14|) becomes 

{u' = F'{u-s)v 

= \(^-eF"{u-s)v^ -v + F{u-s)) . ^' ' 

The singular point u = ^2^, for (|4.14l) corresponds to a pair of equilibrium points (-^2^' 2) 
for dnni), with 



. _ -l±v/l + 4£ F"{zl^-s)F{zl^-7) 
Moreover, from (i) we know that 



U' = Z' - EZ'" = .. . < 0. 



£F'{z-ez";s) 

Hence the orbit lies in the lower halfplane v < Moreover, u — > Z2A: as r/ ^ — cxd, in 
view of the requirements (|4.9() in the original variables and a straightforward analysis of 
the dynamical system (|4.15j) ensures the {z2^-,v\) is a saddle point (see also Section |21 for 
further details on a very close dynamical system). Thus there exists an unique orbit for 
which u{—co) = Z2f^, namely the unique orbit ending at the saddle point. In particular, 
the value v converges toward the finite value < 0, as r] ^ —00. Hence, coming back 
to the original independent variable ^, u' is bounded in 0, that is (|4.13j) . Moreover, as 
already pointed out before, that property implies z"{0) = 0. 

We are now left with the proof of (iii), namely z(^) < '^2k—i' ^'^^ Lemma 12.41 
we shall prove that property by comparing in the phase space {z,z'), the graph (p{z) of 
the trajectory of z{^) with the one ipiz) = F{z;s) of z' = F{z;s) close to the point 
(4„ F(4,; s)). Since V'(4fc) = F'i^k^ = and 

z" 

^\4k) = -7 =0, 
^ 5=0 
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we need to evaluate the second derivative at this point. Clearly, ip"{z2j^) = F"{z2f,;s), 
while 

d (z"\ , 

z'"(0)z'(0) - z"(0)2 z'"(0) 



92 [z^k) = -r- 



(z'(0))3 F(4,;. 



i2- 



?=0 

Moreover, from the above calculations, we know that 
z'(0)-n'(0) _ F{zl^-s)-vl 



z'"{0) 



e e 
Therefore, z"'{^) is the biggest solution of the equation 

Q{\) := (eA - F{z*2,; s)f - = 0. 

Since 

QiFiz*,; sfF"{z*,; s)) = eF[zl^-sfF"{zl^; s)ieFiz*,; s)F"iz*,; s) - 2), 

and F{z2i.; s)F"{z2j^; s) > 0, there exists an eq > such that, for any e < Eq, Q(F(2|^; s)) < 

0, that is 

F{z*2,;sfF"{z*2,;s)<z"'{0), 

which implies 

> ^"iz*2k) > r{zk)- 

Therefore, the graph {p{z) leaves the point Zg^, above tpiz), while going to the left. Hence, 
we can argue as in Lemma 12.41 to conclude z{^) < 22fc-i proof is complete. □ 

In view of the above result, as in Theorem 12.51 we can glue together the profiles zq, 
Zk,h Zk,r, k = 1, ... ,n — 1, Zn to obtain the desired radiating profile, joining n_ = z^ and 
u+ = z+. 

Proposition 4.2. Under conditions \4-4\ ) '^'^^ ( |-<^.5| ), there exists an e > such that, for 
any e < e, there exists a (unique up to space translations) profile z with z(iboo) = z± 
and a speed s such that the function z{x — st) is solution of 6p . This solution is 
away from the n points Z2j^_-^, k = 1, . . . ,n, where z" has at most a jump discontinuity. 

Moreover, there exists a (unique up to space translations) profile u with u(±oo) = u± 
and a speed s (given by the Rankine-Hugoniot condition) such that the function u{x — st) 
is solution of \4-'^ - This profile is continuous away from the n points Z2f._i, k = 1, . . . ,n, 
where it has at most a jump discontinuity. At these points, the Rankine-Hugoniot and the 
admissibility conditions of the scalar conservation law ut + f{u)x = are satisfied. 

Proof. First of all, from Proposition 14.11 we can define profiles zi,. . . ,Zn-i, gluying 
together the profiles Zk,i, z^^r, after an appropriate space translation. Hence, we end up 
with n + 1 profiles zq, . . . , z^, all of them decreasing. 

Moreover, for e < e = mm{ek}k=i,...,n-i, there exist points in the {z,z') plane where 
the graphs of two consecutive can be glued together in that point, in order to give a 
solution of 1)4. 6|) . Moreover, due to the monotonicity of these graphs, again consequence 
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of Proposition 14.11 this intersection is indeed unique. Hence with appropriate space trans- 
lations, we can find a points ^fc, A; = 1, . . . , n — 1, such that Zk{S,k) = ^k+iiCk) = Zk and 
z'kih) = z'k+ii^k) = Zk, k = l,...,n-l. Now we define 

{z„(0 C G (-oo,|'„_i] 
ZkiO, C^[^k,^k-i], k = 2,...,n-l, 
zoiO eG[eo,+oo). 

Thus, this profile defines a solution of H4.6() and it verifies z{±oo) = z±. Moreover, z{x) 
is away from the points ffc, where it verifies ez"{^j. — 0) = = Zfc — /i2fc+2(-Zfc) 

and £z"{ik + 0) = £z'l{ik) = Zk - h2k+i{zk)- 

The regularity oi u = u{x — st) is a direct consequence of the first part of the theorem 
and of the relation u = z — ez" . Moreover, in the case of a discontinuity in ti, namely 
"wfe — 0) / + 0), — 0),n(^fc + 0);s) verifies the Rankine-Hugoniot condition 

for the strictly convex conservation law ut + f(u)x = 0. Indeed, — 0) = /i2fc+2(-2fc), 
u{x + 0) = /i2fe+i(%) and a direct calculation shows 

f{h2k+iizk)) - f{h2k+2{zk)) 

h2k+l{Zk) — h2k+2{Zk) 

In addition, u{^k - 0) = /i2fc+2(5fc) G -^2fc+2, "ufe + 0) = h2k+i{zk) G -^'2^+1 and F{u{^k - 
0); s) = F[u{^k + 0); s) = Since F[-] s) has a local minimum and no local maximum in 
(n(ffc+0),M(^fc-0)), we conclude F(n; s) < F(?x(4±0); s), for any e (u(4+0), M(Cfe-O)), 
that is the Oleinik condition for the inviscid shock {u{^k ~ 0),M(Cfc + 0); s). □ 

Thanks to the above theorem, we know there exist radiative shocks, if Oleinik coindition 
(|4.4|) and the non degeneracy 1)4. 5() are satisfied. Hence, we can pass now to the proof of 
qualitative properties of that solutions. In particular, for a profile u that belongs to BV , 
we shall prove monotonicity and uniqueness in the class modulo L} . It is worth to observe 
that these properties do not require assumption 1)4. 5|) . because they are based on the 
contraction properties of the model under consideration. Actually, they will be used in 
the existence result of Theorem 14.41 in the case of general smooth flux functions which may 
violate 

Theorem 4.3. Let {u,q) be a radiative shock solution of i4-l\ )- -(f ^ G BV, then u is 
monotone. Moreover, given two BV radiative shocks (u, gi) and {v,q2), such that u — v£ 
L^, then u and v are equal up to a space translation. 

Proof. As pointed out before by several authors, |71 IHl El 1131 E] > it is convenient to 
rewrite the 2x2 system ()4.1() as follows: 

ut + f{u)x = -^{u-K'*u), (4.16) 

where the convolution kernel K'^ is given by 

1 1^1 

Let u and v be radiative profiles of class BV, associated with a shock triplet (u_,u-|_;s). 
Assume that v — u is integrable, a fact that is certainly true if u is in BV and v is a shift 
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ThU of u. Finally, let us denote w := v — u. We integrate the entropy inequality 

\w\t + (ifiv) — f {u))sgn{w)) < - {K^ * w — 'w)sgn{w) 
and obtain, with an argument a la Kruzhkov, 

^"rll'"^lli— / {K'^w)sgn{w) dx — \\w\\i. 

Since w is a function oi x — st, the norm ||w||i is constant and the left-hand side above is 
zero. Therefore, we have 

||w||i < / {K'^ * w)sgn{w) dx. 
Jm 

Since has positive values and unit integral, we deduce that w has a constant sign. 

Let /i be a real number. Since ThV := v{- + h) is also a profile and t^v — v £ the 
above fact may be applied to the pair (TfiV,u). We obtain that for every real number h, 
ThV — u has a constant sign. 

Let us focus on the case v = u. Since the integral of t^u — u equals /i(n+ — n_), we see 
that the sign of r/jU — u is that of h(u-^ — U-). In other words, u is monotonous. 

Going back to the general case, the integral of ThV — u equals h{u+ — n_) plus a constant 
(the integral of v — u). Hence there exists a number h for which this integral equals zero. 
But since t^v — u has a constant sign, this means that t^v — u = almost everywhere. □ 

Now we are ready to remove the non degenerate assumption (|4.5() in the existence of a 
radiative shock. We perform this task by means of an approximationg procedure already 
used in ^01 for discrete shock profiles for conservation laws. 

Theorem 4.4. Under conditions i4-4^ , there exists an e > such that, for any e < e, 
there exists a (unique up to space translations) profile z with z(iboo) = z± such that 
the function z{x — st) is a solution of J^.dp , where the speed s is given by the Rankine- 
Hugoniot condition. This solution is away from the n points -^2^-1' k = 1, . . . .,n, where 
z" has at most a jump discontinuity. 

Moreover, there exist a (unique up to space translations) profile u with n(iboo) = u± and 
a speed s such that the function u{x — st) is solution of J^.^ , where the speed s is given 
by the Rankine-Hugoniot condition. This profile is continuous away from the n points 
■^2k-i' k = ■ ■ 7^, where it has at most a jump discontinuity which verifies the Rankine- 
Hugoniot and the admissibility conditions of the scalar conservation law ut + f{u)x = 0. 

Proof. We start by approximating a flux function / satisfying (|4.4jl with a sequence 
of smooth functions /" satisfying both (|4.4I) and (|4.5|) . Then, given an inviscid shock 
(n_,n+;s), n_ > Proposition 14.21 gives the existence of radiative shocks [u^^z^) for 
any n and for any e < E such that u^{-^oo) = z"(iboo) = u±. This value e is independent 
from n, because /" and (/")" remain bounded, (see the constraints that gives the values Ek 
in Proposition ^3 and hence e in Proposition I4.2|l . We fix these profiles with the condition 

^(n_+^.+ )e K(o+),^."(o-)]. 
Since n" is decreasing, 

||r;,n" - u"||ii = \h{u_-u+)\, 
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which imphes that is equicontinuous in L^. Moreover, < u'"' < that is, u'^ is 
equibounded in L°° and hence in Lj^^. We recall that the solution of — e(2:")" + z" = ti" 
with z"(iboo) = u± is given by 

1 /■+°° ICI 

and it satisfies the same estimates. Thus, passing if necessary to subsequences, we have 
u"^ —>■ u and ^ z in Lj^^ and bounded almost everywhere. Clearly, u is decreasing and 
therefore 

u+ < u(+oo) < 2^^^ + u+) < n(— oo) < n_. 
Then, we pass to the limit in the profile equations 
(F^{u^-s) = {z^y 



where F^{u;s) = f^{u) — f^{u±) — s{u — u±) to conclude that {u,z) defines a radia- 
tive profile for the inviscid shock (n(— oo), m(+oo); s). Hence we are left to the proof of 
u(±oo) = u±. Integrating (|4.17j) i and passing into the limit a.e. we get 

z{0-z{0= t F{u{x)-s)dx. 



H 

Since 2:(±oo) = u(iboo) G M and u G from the above relation and from (|4.4|) we 

conclude 

hoo 



0> / F{u{x); s)dx > —oo. 

J — oo 

Therefore in particular F{u{zizoo); s) = and, using (|4.4|) and 

1 



n(+oo) G 



2 {u- + 



, m(— oo) G 



1/ 



we conclude u{±oo) = u±. □ 

Remark 4.5. Proceeding as in SectionlHl it is possible to prove the radiative shocks of (|4.H) 
increase their regularity as e | also in this general non convex case. Indeed, across the 
values -Z2fc~i) ^ = 1, . . . , n, where F{-; s) has local minima, we are in the same situation 
of the one discussed in the convex case and the regularity increases as e | 0. On the 
contrary, across the values Z2i^, k = l,...,n, where F{-;s) has local maxima, we can 
repeat the arguments of Proposition 13.41 and we obtain that the trajectory in the saddle 
point corresponding to the value u = Zgfc is never tangent to the eigenvector (0, 1)* of the 
negative eigenvalue and therefore the profile is regular. 

5. Reduction from the system case to the scalar case 

Let us consider the following strictly hyperbolic-elliptic coupled system 

(ut + f{u)^ + Lq^ = ^^^^ 
I -Qxx + Rq + G ■ = 0, 
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where x S M, t > 0, n G M", q is scalar and R > 0, G, L £ M" are constants. In this section 
we shall prove that the existence of travelling wave solutions of (|5.1|) , for sufficiently small 
shocks, reduces to the study of a scalar model, and therefore we shall obtain their existence 
(and regularity) as corollary of the previous sections. 

According to Definition 11.21 let us consider a radiative shock wave solutions for ()5.1() . 
namely a solution of the form {u, q) = {u{x — st), q{x — st)), with s verifying the Rankine- 
Hugoniot condition 

u(±cx))=u±, s(it+ - n_) = /(n+) - 

As usual, we denote with Ai(u) < • • • < \n{u) the n real eigenvalues of the matrix Vf{u). 
Neglecting the higher order term —qxx in (|5.1|) o (see, for instance jZIlE])) we end up with 
the reduced system 

ut + fiu)x = R~^G ■ UxxL = R~^L Guxx 

and the A:th field has a diffusive dynamics near u±, provided 

4('u±) • L ® Gn{u±) = G ■ rk{u±)ek{u±) ■ L > 0, (5.2) 

where, as usual, r'fc(u) and ik{u) denote the kth right and left eigenvector of V/(n), 
normalized such that £i{u) ■ rj{u) = 6ij. 

Once again, we introduce the variable z as the opposite of the antiderivative of q, that 
is —Zx ■= q and therefore ^;(±oo) = z± = G ■ u±. Hence, we proceed as in the previous 
sections to conclude the the system for u = u{x — st) and z = z{x — st) is given by 



/(n) - /(u±) - s{u - u±) 
z" = G- u. 



(5.3) 



Given the vector L, let P : M" — > and Q : M" — > M be linear applications such that 

kerP = span{L}, QL = 1. 
Then system (|5.c{p becomes 

>(/(n)-/(n±)-s(n-7x±)) =0 

z' = Q{f{u)-f{u±)-s{u-u±)) (5.4) 
G-u = Rz- z". 

Hence, the existence of our profile reduces to a scalar case, provided the constraints H5.4() i . 
together with condition (|5.2j) . ensure 

Q (/(n) - /(n±) - s{u - u±)) = F{G -u-s), (5.5) 

in a neighborhood of {u±; \k{u±)), that is for |n+ — n-_| sufficiently small. This property 
is proved in the next lemma. 

Lemma 5.1. Let f be and assume condition \5.^^) holds. Then, there exists a neigh- 
borhood U X I of {u±; Xk{u±)) such that for any {u; s) £ U x I, 

P{f{u)- f{u±)-s{u-u±)) = ^ u = ^G-u;s), 

with <I> G C"'^(IR X R;Z//), not depending on the choice of P such that kerP = span{L}. 
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Proof. Assume, without loss of generality, (u-i-; Xk{u±)) = and ||G|| = 1. Let Vi, . . . , Vn-i 
be a basis for G-^, so that {Vi, . . . , Vn-i, G} is a basis for R". For u = ciVi+- ■ •+c„_il4i_i + 
aG, let us denote with (/)(ci, . . . , c„_i, a; s) the function P {f{u) — f{u±) — s{u — u±)). 
Then (/)(0, . . . , 0; 0) =0. Moreover, for j E {1, . . . , n — 1}, 

^(o) = PV/(o)y,. 

Then the conclusion follows from the Implicit Function Theorem, provided the vectors 
PV f{0)Vi, . . . ,PV/(0)14_i, j = 1, ... ,n — 1 are linearly independent, because, in that 
case, the relation 0(ci, . . . , Cj_i, a; s) = can be locally written as a function of (a = 
G-U;s). 

Let ai, . . . , ttn-i G M be such that 

aiPVf {0)Vi + ■■■ + a„_iPV/(0)K-i = 0. 

Then 

PV/(0) (aiVi + ■■■ + an-iVn-i) = 
and, by definition of P, there exist /? G M such that 

V/(0) (ai ^1 + • • • + a„__i Vn-i) = PL. 
Applying £fc(0) at the left hand side, we have 

= 4(0)V/(0) (ail/i + • • • + a„_iK-i) = /34(0) • L. 

Hence, condition ()5.2() implies /3 = 0, namely aiVi + • • • + a„_iF„_i € ker V/(0), which 
means 

aiVi H h a„_iK~i = ITkiO), 

for some 7 G R. Since Vi, . . . , Vn~i is a basis for G"*-, using again H5.2|) we conclude 7 = 
and therefore oi = • • • = a„_i = 0. This means that P'Vf{0)Vj, j = 1, ... ,n — 1 are 
linearly independent and Implicit Function Theorem cann be applied. 

Finally, let us note that since kerP = span{L}, P{f{u) — f{u±) — s{u — u±)) = is 
equivalent to f{u) — f{u±) — s{u — u±) £ span {L}, hence the function $ does not depend 
on the specific choice of P. 

The proof is complete. □ 

At this point, we assume the A:th field is genuinely nonlinear, that is 

VAfe(n) • rfc(n) / (5.6) 
and let (u_,u+;s) be a k-Lax radiating shock for ()5.1() . that is 

Xkiu+) < s < Xk{u-), Xk-i{u^) < s < Xk+i{u+). 

We shall prove that in this framework, in the reduced scalar dynamics yielding the exis- 
tence of our profile, that is H5.4() 9. the flux function F{-; s) does not change convexity. On 
the other hand, in the general case and assuming the states verify the Liu E-condition, 
the reduction will lead to the non convex model treated in Section]^ 
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Proposition 5.2. Let us assume condition H5.^) and H5.6^) hold. Then the function 
F{-; s) : M ^ M defined in H5.5}) is either strictly convex or strictly concave in a neighbor- 
hood of z± = G ■ u±. 

Proof. Since (n_,n+;s) form a /c-Lax shock for (|5.1jl . it suffices to prove 



dw"^ 

for sufficiently small shocks. 

Since F{u; s) := /(n) — f{u±) — s{u — u±) = 0, there hold P F{u; s) = and Q F{u; s) = 
F{G ■ u; s) for any u gU. Hence 

F{u;s) = F{G-u;s)L. 

Differentiating with respect to n, one has 

Vf{u)-sI=-^F{G-u;s)L®G (5.7) 
Applying lk{u) and r^iu) respectively to the left and to the right of H5.7() . we get 

Afc(u) - s = -!^F{G-u;s)ik{u)-L®Grk{w) (5.8) 
dw 

Choosing u = u± and assuming |u+ — u_| small enough, the main assumption ^b.2\ implies 
d 



dw 



F{G-u;s) =o(l) as |M_-ti+|->0. (5.9) 



u=u± 



Differentiating (|5.8() with respect to u in the direction of r^iu), we obtain 

VAfc(n) • rfc(u) =^F{G ■ n; s)V(4(n) • L ® Grfc(n)) • rkiu) 
dw 

+ -^F{G-u-s){h{u)-L) (GTfc(n))'. 
Evaluating this relation aX u = u± and taking in account H5.9|) . we obtain 

^MG.u^-s)= ^^^("^)-^^^"^^ +o(l) as 
(4(«)-L)(G.r,(n))^ 

The conclusion follows from GNL condition 1)5.6(1 and assumption ()5.2|) . □ 

The previous results show that the existence of a A;-radiative shock for ()5.ip reduces to 
the study of a scalar model of the form 

{wt + f{w)^+q^ = Q ^^^^^ 
\-qxx + Rq + Wx = 0, 

provided (|5.2() is satisfied for the original shock (ti_,u+;s). In addition, if (|5.6() is also 
verified, we can assume, without loss of generality, that the fiux in (|5.10() is strictly convex 
(see Proposition I5.2() . In this framework, the reduction of (|5.1() to (|5.1UI) is given by (|5.5|) 
and Lemma l5. 11 namely 

w = G-u, F{G • n; s) = Q{f{u) - f{u±) - s{u - u±)), 
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with Q (and P below) as before and, adding if necessary a hnear function to the flux, 

F{w] s) = f{w) - f{w±) - s{w - w±). 

For the sake of clarity, we start by considering the GNL case and we postpone the general 
case at the end of the section. 

If the flux in ()5.1U|) is strictly convex, Theorem l2 . 51 guar antees the existence of a radiative 
shock for that model, with at most a jump discontinuity, which is indeed an admissible 
shock for the inviscid related conservation law. To conclude with the results stated in 
Theorem 11.51 we shall analyze that discontinuity in the corresponding radiative shock of 
the original vectorial case 1)5. If) , showing it forms an admissible radiative shock for that 
system. 

Proof of Theorem \l.f\ Let (n_,n+;s) be an admissible fc-shock for (|5.1|) and, let us con- 
sider = G ■ u±. Let W = ^{U; s) be the neighborhoods of u± and w± given by 
Lemma l5. II and Proposition 15.21 and such that 

4(ti) • Grfc('u) > (5.11) 

for any u £ U. We start by proving that {w-, w+; s) is an admissible shock for the reduced 
scalar conservation law (with strictly convex flux) 

wt + fiw)x = 0. (5.12) 
Indeed, from 

f(.u+) - f{u^) - s{u+ - U-) = 0, 
we obtain in particular 

F(u;+; s) = Q{f{u^) - f{u.) - s{u+ - n_)) = 0, 

that is, the Rankine-Hugoniot conditions of 1)5. 12() for the shock {w-,Wj^] s). 
Proceeding as in the proof of Proposition 15.21 we obtain 1)5. 8p . that is 

\k{u) - s = {f'{G ■ u) - s)£kiu) ■ L (S)GrkH- (5-13) 

Therefore the sign of f'{G -u) — s is given by the sign of Afc(u) — s for any u £U, in view of 
1)5. 11(1 . Thus, using 1)5. 13() for u = u± and taking into account the admissibility condition 
of u±, AA;(n+) < s < Afc(n_), we obtain 

f\w.)>s>f'{w+), 

that is, the discontinuity {w-,w+] s) is admissible for ()5.12() . 

At this point, let w be the (unique up to space shifts) radiative profile of the reduced 
model (|5.1Uj) given by Theorem 12.51 Then, the above results guarantee the existence of 
the (unique up to space shift) radiative profile {u,q) for ()5.1|) . if |n_ — is sufficiently 
small. Therefore we only have to prove that, if w is discontinuous, the corresponding 
discontinuity in u defines an admissible /c-shock wave for the inviscid hyperbolic system 
of conservation laws ut + f{u)x = 0. We shall perform this task as before, connecting the 
properties of the shock for the reduced system with the ones of the shock of the original 
systems. 
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Denoting with {wi, Wr',s) the discontinuity of the radiating shock for the scalar reduced 
model, we have wi^r S W and therefore there exist unique ui^r £ ^ solutions of 

\p(f{u)-fiu±)-siu-u±)) = 0, ^ ■ ' 

which are defined through the function <I> constructed in Lemma l5. II 

As before, let us start by showing the Rankine-Hugoniot condition for (u/, u^; s), namely 

fiui) - f{Ur) = S{ui - Ur). 

Clearly, the above relation is equivalent to 

F{ui;s) = F{ur;s). (5.15) 
In order to prove (|5.15j) . we only have to prove 

QF{ui-s) = QF{ur;s), (5.16) 

because PF{ui; s) = = PF{ur\ s) is given by 1)5. 14() '?. Moreover, QF{u; s) = F{G ■ u = 
w; s) for any u €U, namely, for u = ^{w; s) and therefore ()5.16|) is precisely the Rankine- 
Hugoniot condition 

f{wi) - f{Wr) = s{wi - Wr) 

for the reduced scalar model (|5.1U|) . 

We turn now to the proof of the Lax conditions for the A:-shock {ui,Ur] s) of ut + f{u)x = 
0, namely 

\k{Ur) < s < Xk{ui), (5.17) 
Xk-l{ui) < S < Xk+l{Ur). (5.18) 

Since we are dealing with weak shocks and the system is assumed to be strictly hyperbolic, 
(|5.18() follows from (|5.17|) . Finally, using once again 1)5. 8() . this time for u = u^r^ and taking 
into account the admissibility condition of wi^r = G ■ ui^^i that is 

f'{wi) > S> f{Wr), 

we obtain (j^T7|) . □ 

Proof of Theorem ] l.<^ The proof is an immediate consequence of Proposition mi and The- 
orem 11.61 □ 

Let us pass now to the general case, namely when condition ()5.6() is violated. As we 
have shown in Section^ the existence of radiative shock without convexity assumptions 
is guaranteed only when the radiative effect is sufficiently dissipative. Hence, we shall 
consider the following system 

(ut + f{u)x + Lqx = 
\-eqxx + Rq + G ■Ux = 0, 
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with < e ^ 1. It is worth to observe that the smallness in e is needed only for the 
existence of the profile for the reduced scalar model 

U + Kw)^+q. = ^^^^^^ 
[-eq^x + Rq + = 

and it does not play any role in the connections between the admissibility conditions for 
the jumps of that model and the ones of ()5.19|) . 

Theorem 5.3. Let {u-,u^; s) be an admissible k -shock for M.^) and assume \5.i3fl holds. 
If\u^ —u^\ and e are sufficiently small, then there exists a (unique up to shift) admissible 
radiative shock wave {u = u{x — st),q = q{x — st)) of \5.iy\) such that (n(±oo), g(±oo)) = 
(n±,0). 

Proof. The proof of this theorem follows the same lines of the one of Theorem 11.61 In 
particular, the existence of a radiative shock for (|5.19j) comes from Theorem 14.41 and the 
reduction to 1)5.20(1 . Clearly, the analysis of Rankine-Hugonoit conditions is made as 
before, because it is independent from convexity assumptions. Theorefore we are left to 
the proof that, given the neighborhoods U and W as before, a shock («;,itr;s), ui^r G U, 
verifies the Liu E-condition for (|1.2|) if and only if the corresponding shock (wi,Wr',s), 
wi^r G W, verifies the Oleinik condition for 

wt + fiw), = 0. (5.21) 

Let s) be a A:-shock of (|5.1() which verifies the (strict) Liu E-condition: denoting 

with Uki^r) the k-th shock curve, U- = tffc(O), ti+ = Uk{f), then 

s = Sfc(f) < Skir) (5.22) 

for any r between and f . As before, consider w± = G ■ u± and assume W- > w+. By 
definition of shock curve, we have 

= f{Uk{T)) - f{u_) - Sk{T){uk{T) - U_) 

= F{uk{T); s) + {s- Sk{T)){uk{T) - n_) 

= F{G ■ Uk{T)]s)L + (S - Sk{T)){Uk{T) - U-) 

which implies 

F{G ■ Ukir); s)L = (s,,(t) - s)(n,,(T) - u.). (5.23) 

Since for small shocks, that is |r| small, Ukir) — U- = Tr^iu^) + o(r), form ()5.23|) we 
conclude 

F{G ■ Uk{T); s)L = (r r,,(n_) + o(T))(s,,(r) - s) 

and, multiplying that relation on the left for ik(u-) 

F{G ■ Uk{T); s)ik{u^) ■L = {t + o(T))(sfc(r) - s). (5.24) 

Moreover, we differentiate Wkir) := G ■ Uk{T) with respect to r to conclude Wkir) = 
G-Uk{T) = G ■ rk{u-) + o{t) , and, since we are dealing with small shocks and G-r^iu ) ^ 
for (|5.2() . it follows that wi^^t) is decreasing for r between and f. Hence, if G-rk{u-) > 
(resp. G ■ ri^{u^) < 0), then f < (resp. r > 0) and ik{u-) ■ L > (resp. £fc(ti-) • L < 0) 
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using again (|5.2|) . Therefore, for |t| small and between and r, r + o(r) < (resp. 
r + o(r) > 0) if 4(n_) • L > (resp. 4(n_) • L < 0) and therefore from (|^^ and (|^?^ 
we conclude F{G ■ Ufc(r);s) < for any r between and f, that is, the strict Oleinik 
condition for the shock {w-,w^; s) of H5.2H) . 

Let us now consider a radiative profile for ()5.2U|) which has a discontinuity {wi,Wr;s) 
admissible for (|5.21j) . As in the proof of Theorem ll.6l for the GNL case, there exist unique 
ui^r in the neighborhood under consideration such that the discontinuity {ui,Ur; s) verifies 
the Rankine-Hugoniot conditions of (|1.2() and in particular it belongs to a fc-shock curve 
Uk{T) with Uyfc(O) = ui, Uk{Tr) = Ur- The relation 

f{uk{T)) - f{ui) - Sfc(r)(ufe(r) - m) = 
implies this time 

{F{G ■ Ufc(r); s) - F{G ■ u«,r(T); s)L = (sfc(r) - s)(ufe(r) - U-), 
which becomes for |r| small 

{F{G ■ Ukir); s) - F{G ■ nz,,(r); s))ik{ui) ■ L = (sfc(r) - s)(r + o(r)). 

Moreover, using (jS-llf) this time for ui, we argue as before to obtain the (strict) Liu 
E-condition 

S = Sk{Tr) < Skir) 

for any r between and from the (strict) Oleinik condition 

F{G-Uk{r);s) <F{G-ui,r{Ty,s) 
for any Wkir) = G ■ Ukir) £ {wr-,wi) and the proof is complete. □ 
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